A nonlocal wave-wave interaction among Alfven waves in an intermediate-/? plasma 
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A nonlocal coupling mechanism to directly transfer the energy from large-scale Magnetohydro- 
dynamic(MHD) Alfven waves to small-scale kinetic Alfven waves is presented. It is shown that 
the interaction between a MHD Alfven wave and a reversely propagating kinetic Alfven wave can 
generate another kinetic Alfven wave, and this interaction exists in the plasmas where the thermal 
to magnetic pressure ratio is larger than the electron to ion mass ratio. The proposed nonlocal in- 
teraction may have a potential application to account for the observed electron scale kinetic Alfven 
waves in the solar wind and solar corona plasmas. 

PACS numbers: 



I. INTRODUCTION 

The large-scale magnetodydrodynamic(MHD) Alfven 
wave 1 -, which is generated by the magnetic stress tensor, 
is a nondispersive low frequency wave mode and can be 
directly derived from MHD equations. The dispersive 
Alfven wave arises due to the finite-frequency and small- 
scale modifications 2 . The dispersive Alfven wave is spe- 
cially called as the kinetic Alfven wave (KAW) when its 
perpendicular wavelength is close to the order of the ion 
gyroradius 3 -, the ion acoustic gyroradius, or the electron 
inertial lengthl The KAW can accelerate and heat elec- 
trons by its parallel electric field 5-7 and interact with ions 
through its perpendicular electric field^— . The proper- 
ties of the KAW have been demonstrated by many exper- 
imental investigations in spaced— and laborotarjsi^— 
plasmas. In particular, the KAW can play an important 
role in the acceleration of energetic electrons in Earth's 
aurora—"— and solar flares 2 ^— as well as the heating of 
solar coronal plasmas 2 ^—. 

Extensive works have been done to discuss the local 
interaction process of three Alfven waves in which the 
magnitudes of three wavenumbers are comparable^ 2 . - — . 
A large amplitude dispersive Alfven wave can bring about 
the local nonlinear decay among themselves^ 4 . - —. Local 
interaction between two counterpropagating Alfven wave 
packets can lead the wave energy to cascade from the en- 
ergy injection region to the energy dissipation region in 
the MHD turbulenco 32 ' 33 ' 38 ri£. Through the local cas- 
cade, the KAW can be generated at scales of the order of 
the ion inertial length or the ion gyroradius in the solar 
wind turbulence 41 - - — , however, it cannot reach electron 
scales (the electron inertial length or the electron gyrora- 
dius) due to the large electron Landau damping in these 
small scales— 

Recently, some works have found that the nonlin- 
ear coupling of Alfven waves with different scales can 



occu r 46 ' 47 . Voitenko and Goossens^ 6 - presented a nonlo- 
cal interaction among three Alfven waves and showed 
that a large-scale MHD Alfven wave can decay para- 
metrically into two small-scale KAWs. Shukla and 
Stcnflc 47 studied the three-wave interaction involving 
two KAWs and one field-aligned Alfven wave with the 
finite-frequency modification, and they also showed two 
small-scale KAWs can be nonlocally excited by this large- 
scale Alfven wave. Unlike the local energy cascade in the 
MHD turbulence, the wave energy can be directly trans- 
ferred from large-scale Alfven waves to small-scale KAWs 
in these nonlocal interaction processes. 

In this paper, we investigate the nonlocal interaction 
among one MHD Alfven wave and two KAWs, and we 
shall show that the nonlocal coupling, a MHD Alfven 
wave + KAW —> KAW, may play an important role in 
generating the KAWs with electron scale a 43 ' 44 . The re- 
minder of this paper is organized in the following fashion. 
The qualitative and quantitative analyses are given in 
section 2 and section 3, respectively. Section 4 presents 
an application in the solar corona. The discussion is set 
in section 5 and the summary is contained in section 6. 



II. QUALITATIVE ANALYSIS 

Three waves in the nonlinear coupling process (a MHD 
Alfven wave + KAW 1 — > KAW 2) must satisfy the reso- 
nant relation, which describes the phase relation of these 
three waves, 



- Ul\ = UJ-2 

ki = k 2 



(1) 



* Corresponding author. Electronic address: djwuapmo.ac.cn 



where uj s and k s are the frequency and wave vector 
of the large-scale MHD Alfven wave, respectively; W1.2 
and ki^ are frequencies and wave vectors of the two 
KAWs, respectively. In this study, we consider the gen- 
eral oblique propagating MHD Alfven wave, where wave 
vector k s = k s ^ + fc S2 z, and k s j_ and k sz z are the wave 
vectors perpendicular and parallel to the background 
magnetic field -Bqz, respectively. 
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The linear dispersion relation of the MHD Alfven wave 
is uj s = Va^-szi where Va is the Alfven velocity. For the 
two KAWs in ([T]), their linear dispersion relations are 

given as uip = Vkfci,2»i^i,2 with Kx,2 = yl + P 2 ^i,2J_ 
for the intermediate-beta plasmas (m e /mj <C /3 <C 1) 
and if ! 2 = 1/yl + ^e^i 2_l ^ or ^ ne low-beta plasmas 

(/3 -C m e /m,i)2^, where fci :2z and fei : 2_L are the parallel 
and perpendicular wavenumbers of the two KAWs, re- 
spectively, p 2 — p 2 + p 2 {pi is the ion gyroradius and p s 
is the ion acoustic gyroradius), A e is the electron inertial 
length, f3 is the ratio of the plasma thermal to magnetic 
pressures, and m e /mi is the electron to ion mass ratio. 

Substituting the linear dispersion relations of the MHD 
Alfven wave and two KAWs into the frequency relation in 
(fl]) and then combining the z-direction wave-vector rela- 
tion, we can obtain a restricted relation for two wavenum- 
bers of the two KAWs, 
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FIG. 1: The parallelogram in the (cj, k z ) plane for the three 
Alfven waves interaction, a MHD Alfven wave + KAW 1 — > 
KAW 2. 



k\ z {Ki - si) = k 2z (K 2 - s 2 ), 
or be written as 

(i^i-si)(K 2 -s 2 ) >0. 



(2) 



(3) 



where the MHD Alfven wave is assumed along the back- 
ground magnetic field. The subscribes Si ; 2 denote the 
directions of the two KAWs, for example, s\ t2 — 1 stand 
for two waves having the same directions as the MHD 
Alfven wave and Sx,2 = — 1 are two reversely-propagating 
waves. 

There exist three kinds of wave-wave interaction, 
Si t 2 = ±1 and s\ — — s 2 = — 1, for the coupling of 
one MHD Alfven wave and two KAWs. In the first 
two interaction cases (si.2 = ±1), the frequency of the 
KAW 1 is required much larger than the MHD Alfven 
wave frequency, and we do not discuss these two cases 
in this study. The third case (si = — s 2 = —1), 
as shown in Fig. 1, describes the parallel-propagating 
KAW 2 generated in the nonlinear interaction between 
the parallel-propagating MHD Alfven wave and the 
reversely-propagating KAW 1. From the restricted rela- 
tion Eq. (J3J, we can see that the third interaction takes 
place under the condition K 2 > 1, which means this in- 
teraction only exists in the plasmas with j3 > m e /mj. 
Further, we can obtain u>i — u> s (K 2 — \) / (K 2 / K\ + 1) 
and uj 2 = w s (K\ + 1)/(K 1 /K 2 + 1) from Eqs. Q and 
©■ 



III. TWO-FLUID MODEL 

Let us consider the homogeneous and collisionless plas- 
mas with two species of particles, electrons and ions. 
We restrict the plasmas in the intermediate-beta range 
(m e /mi <C |3 < 1), where the ion finite Larmor radius 
and ion polarization contributions are both important for 



the KAWs^. The momentum equation is, 



dtv r 



T, 



-^-(E + Vj xB ) 

m,- 



3 
«7 



Vj x B, 



(4) 



where Vj, mj, q.j, nj, and Tj denote the particle veloc- 
ity, mass, charge, number density, and temperature for 
species j, respectively; E and B are the wave electric and 
magnetic field perturbations, respectively. Note that the 
isothermal assumption has been used in (4). 

The variables can be written as, v 3 = v JS + Vji + Vj 2 , 
rij = n Q + riji + nj2, E = E s + Ei + E 2 and B = 
B s + Bi + B2 (where no is the background number den- 
sity) . The large-scale MHD Alfven wave only has the per- 
pendicular velocity and electromagnetic perturbations: 
v s = v s j_, E s = E s _l and B s = B s ^. By substituting 
these expressions into the momentum Eq. we can 
obtain the fluid velocities of the KAW 2, 



1 v Tirj n a 

v i2 _l ^ ^-E 2 ± x z —V± x z - 

iso Wei no 



1 



BqLO c 



-9*E 



n i2 
n 



■\y ia ± ■ Vj_(vii£ + v ilD )] x z, (5) 



1 %r7 n e2 

v e2 _L -^-E 2i x z — V_l x z 

+ 4~( V elz X B s-L) X Z. 
t>0 



(6) 



and 



e „ 2 q n e2 
■fUeSz — &2z - V T d z 

m e no 



dm 



[(VelE + VelD) X B s j_ + V es j_ X Bij_] • Z. 

m e 



(7) 
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where uj C i — eBo/rn-i and uj ce — —eBo/m e are the 
ion and electron cyclotron frequencies, respectively; 
vn = {Ti/mi) 1 / 2 and vt £ = (Te/m,,) 1 / 2 are the ion 
and electron thermal velocities, respectively; Vhe — 
(1/S )Ei_l x z and v aD = -(v^Ju ci )V ±(n a /n ) x 
z denote the ion electric field drift velocity and the 
ion diamagnetic drift velocity of the KAW 1, respec- 
tively; VeU, VelE = (1/-B )Ei_l x z and \ elD = 
— (vT e /u ce )V±(n e i/n ) x z denote the electron parallel 
velocity, the electric field drift velocity, and the electron 
diamagnetic drift velocity of the KAW 1, respectively. 
Here, the nonlinear terms due to the thermal pressure 
gradient are neglected because they are smaller than the 
nonlinear effects in ((5|) - ((TJ). This study is limited to 
the weak nonlinear system, where a Fourier analysis can 
be used. That is, the wave amplitudes are restricted as 
small amplitudes, and the ratio of energy in the oscilla- 
tions to the total energy of the plasma is a small param- 
eter (e.gi^S). The waves in this system are assumed to 
be three monochromatic waves, and their amplitudes can 
be expressed as Ai = A kl e~ l " )lt+l ' kl ' r , where / = (l,2,s). 

The ion and electron number density perturbations in 
equations © - can be obtained by using the ion 
continuity equation and the quasi-neutrality condition 
(n 2i = n 2e = n 2 ), 



(1 + ptkl ± )n 2 /n 



-k2J_ • E 2 j_. 



(8) 



The wave electromagnetic fields E 2 and B 2 in J5]) - J8]) 
can be conveniently represented by a scalar potential <p 2 
and a z-direction vector potential A 2z z, 

E 2 = -ik 2 2 + iuJ 2 {A 2z z), B 2 = ik 2± x (A 2z z). (9) 

By using (J5j) - (JSJ) , the current density J 2 = eng (v.i 2 ± — 
v e 2_L — v 2ez z) is written as, 



1 



;2 ± n oe 
-uj 2 k 2± (p 2 k 2 _i_- 



{[v is j_ • ikix(Vii B + V s ;i_d) + — V elz X B sJ _] x z}, 

rrii 



(10) 



and, 



n e 2 l + p 2 k 2 ± 
0O 2 J 2 z — — — 2T2 - s 2 k 2z (p 2 



n e 



m e 1 + pjk£ 



oj 2 A 



2z 



n e 



[(VelE + V g id) X B s± + V es J_ X BiJ_] • Z, (11) 



The current density can also be given in terms of A 2z 
by Ampere's law po3 2 = V x B2, 

MoJ2_l = -k 2 ±k 2z A 2zi fi Q J 2z = k 2X A 2z , (12) 

From Eqs. (fit))) and (fT2")l . a nonlinear equation for the 



KAW 2 is given by, 

[uJ 2 2 -V^l + p 2 k 2 2± )k 2 2z ^ 2 
- B (l + pf fc| ± )w 2 fc~_fk 2± 



{[v is ± ■ iki_L(viiB + Vni 



-V e u X B s _l] X z} 



+ Vl{\ + p 2 k 2 21 )is 2 k 2z 

[{VelE + Veil?) X B s± + V esJ _ X BijJ • Z. (13) 

where the first nonlinear term on the right-hand of the 
equation comes from the ion convective motion and per- 
pendicular electron nonlinear Lorentz force, and the sec- 
ond nonlinear term arises due to the parallel electron non- 
linear Lorentz force. Three nonlinear effects (the ion con- 
vective motion, the parallel and perpendicular electron 
nonlinear Lorentz forces) are all important in the nonlo- 
cal coupling of one MHD Alfven wave and two KAWs, 
and these three effects are also contained in some other 
nonlinear models in simplified forms (e.g. , 37 ' 47 ' 49 ). 

The variables of the KAW 1 in the nonlinear Eq. (fT3|) 
can be expressed in terms of the scalar potential 4>\ , 



Viie = v elE = -— kix x z^i, 
i p 2 k 2 \ 

VilD = — ; ~ 2 u2 kl -L X Z ^ 1 ' 



Velz 



B Q 1 + pfkf ± 

i p 2 s k 2 ± 



Bo 1 + p\k{ 

I U2 



2L2 X 
k 2 ± SlCJi 



B lo ci 1 + pfkf ± fci 
. si^i k 1± x z 

Bl± = Vj^i + p^> 



(14) 



For the obliquely propagating MHD Alfven wave, we 
use a scalar potential <p s to express the electric field per- 
turbation E s = —V±4> s and other variables in Eq. (fTB"]) . 



Va 
B 



B, 



Bo 



k s± x z0 s . (15) 



By using linear relations in Eqs. (fT4"l) and (fT5)) , we can 
obtain the nonlinear dispersion relation of the KAW 2, 

±>o i + p { k 1± 

(sis 2 ^ + -^i)(l - si#i)(k s± x z) • \c 1± Ms- (16) 
K 2 k£ ± 



where K\_ 2 = yl + P 2 fci2± * s the definition in the 
intermediate-beta plasmas and the relation k\± ~ k 2 ± 3> 
k s j_ has been used. The nonlinear dispersion relation of 
the KAW 1 can be obtained by exchanging the subscripts 
1 and 2 and replacing <fi s by <fi* s in Eq. (fTB)) . We here re- 
strict the MHD Alfven wave as the finite amplitude wave 
and neglect its nonlinear modification caused by the two 
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KAWs. This finite amplitude assumption can be permit- 
ted for some solar and space plasma environments, such 
as the solar corona and the solar wind. 

Two KAWs frequencies are expressed as u>i — u>i r — iji 
for the KAW 1 and W2 = u 2r + 1J2 for the KAW 2, where 
w r i,2 denote the real parts of two frequencies, 71 is the 
decay rate of the KAW 1 and 72 is the growth rate of the 
KAW 2. To show which parameters affect the nonlocal 
coupling strength among three Alfven waves, we assume 
71 ~ 72 = 7 and 7 <C oj r . Then, 7 can be derived from 
the nonlinear dispersion relations of the two KAWs, 

2 _ V% K 2 k 2 2± 
' 4 K\ 

(siS2^ + Sh 2 ^! - si)(K 2 - S2 )sin 2 ^, (17) 

where B s ± — k s ±(j} s /VA and is the angle between k s j_ 
and kix. From the equation (jTTJ) , the restrict relation 
(A'i — si)(A2 — S2) > can be reobtained. The equation 
(jTTJ) is rewritten as the following form for the interaction 

case s\ = — S2 = — 1, 

a _ V\ K 2 k* ± 
7 4 K x 

K 2 k 2± B 

IV. APPLICATION IN THE SOLAR CORONA 

Here we use some typical parameters in the quiet so- 
lar corona to discuss the nonlocal interaction proposed 
in this study: the number density n = 10 9 cm~ 3 , the 
thermal temperature Tj = T e = 10 6 K, and the mag- 
netic field Bq — 10 G. These values give the plasma 
beta j3 ~ 3.5 x 10~ 2 that is located in the intermediate- 
beta range. The existence of the MHD Alfven wave 
has been demonstrated in the solar corona plasmas by 
many observations^—. The frequencies of these MHD 
Alfven waves are in the range (10 -5 Hz, 10 -1 Hz)£2r— , 
or the high frequency range (0.1 Hz, 2.5 Hz) if these 
waves are generated by the reconnection in the magnetic 
networks^. As an example, we consider two MHD Alfven 
waves with the same frequencies 10 _1 Hz but different 
directions, k s ± = k sz (pk s ± ~ 10 -7 ) and k s ± = 10k sz 
(pk s ± ~ 10~ 6 ). If the two KAWs locate in the frequency 
range (10 -5 Hz, 2.5 Hz), we can use the frequency rela- 
tion u>i = uj s (K2 — I) / (K2/ K\ + 1) to estimate the limi- 
tation of the perpendicular wavenumber for the KAW 1. 
It shows that pki± is required larger than 0.02, mean- 
while, we can also obtain the limitation of pk 2 ±_ because 
of pki ± ~ pk 2 ±. 

Figure 2 gives the relation between the maximal cou- 
pling strength 7 and the perpendicular wavenumber of 
the KAW 2 pk2±, where the solid line denotes the MHD 
Alfven wave propagating case k s ± = k sz and the dashed 
line is for the case k s ± = 10k sz . Figure 2 shows that the 
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FIG. 2: The relation between the maximal coupling strength 
7 and the perpendicular wavenumber of the KAW 2 pk2±- 
The solid line denotes the MHD Alfven wave propagating case 
k s ± = k sz and the dashed line is for the case k s ± = 10k sz . 



coupling strength 7 increases with the increment of the 
pk2± and the 7 is large for the more obliquely propagat- 
ing MHD Alfven wave. 

The new generated KAW can suffer the electron Lan- 
dau damping in the collisionless plasmas^, and the ex- 
pression of the electron Landau damping rate is given 

as 7l/w cj; = A(m e /m i ) 1 / 2 l3l /2 (X i k ± ) 2 X i k z ^; where /3 e 
is the electron plasma beta, Xi is the ion inertial length 
and A is a function of /3 e and T e /Tj. For the low beta 
plasmas with f3 e <C 1, A is roughly 0.4. We can use the 
electron Landau damping rate and the expression (|18[) to 
give some estimations for the threshold amplitude of the 
MHD Alfven wave. 

Figure 3 gives the threshold amplitude of the MHD 
Alfven wave as a function of the pk2±. Figure 3 shows 
that the threshold amplitude is smaller for the more 
obliquely propagating MHD Alfven wave and for the 
KAWs with smaller pk2± . In the solar low corona, the rel- 
ative amplitude of the MHD Alfven wave nearly reaches 
0.2^2, which is larger than all the threshold amplitudes 
shown in Figure 3. A small threshold amplitude, espe- 
cially for the more obliquely propagating MHD Alfven 
wave, implies that the nonlocal interaction mechanism 
may happen in the solar corona. 



V. DISCUSSION 

Our nonlinear wave-wave interaction requires a large- 
scale MHD Alfven wave and a small-scale KAW with 
counterpropagating directions. These two waves can co- 
exist in the solar and space plasmas. The MHD Alfven 
wave is a common wave mode and can be easily gen- 
erated by the shear motion and/or the magnetic recon- 
nection process. The KAW can be produced by many 
mechanisms, such as the phase mixing and the resonant 
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tron gyroradius scale) in the solar win d 43 ' 44 ' 65 . This elec- 
tron scale KAW cannot be explained by the local wave 
cascade mechanism because the large electron Landau 
damping can totally dissipate the KAW energy in that 
small scale region 4 ^. But our results implies, except for 
the local wave-wave interaction, the nonlocal wave-wave 
interaction may also happen in the Alfven wave turbu- 
lence. Furthermore, as shown in figure 2, the nonlocal 
coupling is stronger for the smaller scale KAWs. So the 
observed electron scale KAWs may be produced by this 
nonlocal mechanism. It should be pointed out that the 
nonlocal turbulence of the KAWs still receives little at- 
tention and this problem will be further investigated in 
our future works. 



FIG. 3: The threshold amplitude of the MHD Alfven wave 
changes with the perpendicular wavenumber of the KAW 2 
The lines have the same meaning as in Fig. 2. 

absorption of the MHD Alfven -waves**—, the instabil- 
ity caused by the warm proton beams or by the ion-ion 
streamin g 60 ' 61 , the velocity shear transition^ 2 -, and the 
parametric decay by other wave mode o 63 ' 64 . These wave 
generation mechanisms make it easily to produce the op- 
posite traveling MHD Alfven wave and KAW, thus caus- 
ing the nonlocal interaction. 

The mechanism proposed in this study is also dif- 
ferent from the other two kinds of interactions among 
three Alfven waves, the local interaction and the non- 
local deca y 32 ' 33 ' 38-40 ' 46 ' 47 . The proposed wave-wave in- 
teraction happens for Alfven waves with different scales, 
whereas the local interaction is limited to waves with 
comparable scale s 32 ' 33 ' 38 " — . The nonlocal decay mecha- 
nism also relates to three waves with different scales, how- 
ever, it describes two waves excited by a pump wav o 46 ' 47 . 
Of course, two nonlocal interactions, either through the 
nonlocal decay or through our nonlocal interaction, can 
both directly transport the Alfven wave energy from the 
large-scale region to the small-scale region and then dis- 
sipate the wave energy there. 

Last, let us give a simple discussion for a potential 
application of our nonlocal mechanism. Recent obser- 
vations showed there exists the KAW even down to the 
electron scales (the electron inertial length or the elec- 



VI. CONCLUSION 

This paper discusses the nonlocal coupling of one MHD 
Alfven wave and two KAWs, a MHD Alfven wave + 
KAW 1 — > KAW 2. The qualitative discussion shows 
that this wave-wave interaction works in the plasmas 
with /3 > m e /mi. The frequency relations among the 
two KAWs and the MHD Alfven wave are ui\ = u> s (K2 — 
1)/(K 2 /K 1 + 1) and lo 2 = u) s {K x + 1)/{K X /K 2 + 1). 
The quantitative discussion shows that the coupling 
is stronger for the more obliquely propagating MHD 
Alfven wave and for the KAWs with lager perpendicular 
wavenumbers. We also show that the proposed mecha- 
nism may play an important role in generating electron 
scale KAWs in the solar wind. 
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